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We present protocols for implementation of universal quantum gates on an arbi- 
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. trary superposition of quantum states in a scalable solid-state Ising spin quantum 
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£NJ \ computer. The spin chain is composed of identical spins 1/2 with the Ising inter- 

CN ■ action between the neighboring spins. The selective excitations of the spins are 
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. provided by the gradient of the external magnetic field. The protocols are built 

^ i f , of rectangular radio-frequency pulses. The phase and probability errors caused by 
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' unwanted transitions are minimized and computed numerically. 
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I. INTRODUCTION 

Current proposals on a scalable solid-state spin quantum computer can be roughly divided 
into two main streams. The first one relies on a controllable interaction between the electron 
spins (see, for example, [l], @). In such proposals a point contact gate switches interaction 
between the neighboring qubits. It was shown, in particular, that the switchable Heisenberg 
interaction alone can provide the universal quantum computation. The main advantage 
of the controllable interaction proposals is the high clock speed in the GHz range. The 
main technical challenge in such proposals is the precise control of the interaction. It was 
estimated that the interaction constant must be controllable with the accuracy 10 -4 . Also 
the time of turning on (off) the interaction must be about 10-100 fs. The other challenge is 
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a possible decoherence caused by the point contact. 

The second stream relies on the permanent interaction between spins (see, for example 
@, U). m sucn P ro P osa l s electromagnetic pulses provide implementation of quantum com- 
putation. The main advantage of these proposals is the use of well developed technique of 
electromagnetic pulses, and the absence of the electro-static gates. The main technical chal- 
lenge for these proposals is the creation of a large gradient of the magnetic field. The other 
challenge is the decoherence caused by the sources of the magnetic field. The latest devel- 
opment of the micropattern wires technique provides the magnetic field gradients 10 5 — 10 6 
T/m [||, ||. Even greater gradients ~ 10 7 T/m are expected if the micropattern wires will 
be replaced by the ferromagnetic plates 0. This development provides new ground for the 
electromagnetic pulse proposals. The field gradient 10 6 T/m corresponds to the magnetic 
field difference 0.03 T for the distance 30 nm. In turn, the difference between the electron 
spin resonance (ESR) frequencies Au/2tt for this distance will be 840 MHz. 

For the electromagnetic pulse proposals the most convenient interaction between qubits 
is, definitely, the Ising interaction [§, ||. Normally, the Ising interaction is considered as 
non-typical for solids. However, if the frequency difference between the neighboring spins 
is greater than the spin-spin interaction, the spin-spin interaction (Heisenberg or dipole- 
dipole) becomes effectively the Ising interaction. This effect is well-known in liquid nuclear 
magnetic resonance (NMR) where the scalar coupling, which is equivalent to the Heisenberg 
interaction, becomes the Ising interaction due to the chemical shift of the NMR frequencies 



10| , PHI . The same effect takes place for the heteronuclear dipole-dipole interaction in 



solids. We should note, that the long-range dipole-dipole interaction in a spin chain can be 
effectively suppressed if the angle between the chain and the external magnetic field equals 
to the magic angle fTTJ, |j~2"f . 

We believe that the large field gradients provide the resurgence of the interest to the Ising 
spin quantum computer. If the interaction between the neighboring paramagnetic spins is of 
the order of 10MHz and the ESR frequency difference is Aa;/27r = 840 MHz, the spin-spin 
interaction becomes effectively the Ising interaction. As an example, we consider phosphorus 
impurity donors in silicon. In the first approximation, the exchange interaction constant for 
the donor electrons separated by the distance r is given by [IJ 

- 2 , r ^5/2 



e 2 /r \5/2 
J(r) = 0.8 — - exp(-2r/a ). 
ean Van/ 
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Here e is the dielectric constant of silicon, a® is the effective radius of the impurity atom, and 
e is the electron charge. Taking e = 12, ao = 3 nm, and r = 30 nm, we obtain J = 5 MHz. 

Currently, there are three main approaches for implementation of quantum logic gates 
using the electromagnetic pulses. The first approach is developed in experiments with liquid 
state NMR quantum computer (see, for, example, JTIJ). To provide the conditional logic 
this approach takes advantage of the spin-spin interaction during the time-interval between 
the electromagnetic pulses. Thus, the Rabi frequency Q of the pulses [Q = jBi, where 7 
is the gyromagnetic ratio and B\ is the amplitude of the radio frequency irf) field] is much 
greater than the Ising interaction constant J, but the clock speed of the gate is J. 

The second approach is intended to increase the clock speed by application of more 
powerful strongly modulated pulses which average to zero all undesired interactions between 
spins but drive the desired quantum transitions. The first experiments indicate the almost 
tenfold increase of the clock speed Jl4| . 



The third approach is based on the application of the selective electromagnetic pulses with 
the Rabi frequency Q < J (see, for example, [0, [15], [Tj|). These pulses drive a resonant spin 
depending on the states of its neighbors. The clock speed in this approach is determined by 
the Rabi frequency, fl While the third approach has a clear disadvantage in the clock speed, 
its advantage is the small power of the pulses: for a given angle of the spin rotation the 
power of the pulse is inversely proportional to its duration. The small power of the pulses is 
not important for the room temperature NMR quantum computer but can be important for 
a scalable solid-state quantum computer, which is expected to operate at low temperature. 
Besides, the powerful pulses can contribute to the decoherence rate. 

Based on the above remarks, we believe that the development of the theoretical back- 
ground for the quantum logic gate implementation using the selective electromagnetic pulses 
in the Ising spin quantum computer is an important task. In this paper, we present a scheme 
for implementation of universal quantum gates on an arbitrary superposition of quantum 
states in the Ising spin quantum computer. 

The paper is organized as follows. The general problem to be solved is formulated in 
Sec. II. The Ising spin quantum computer model is described in Sec. III. The probability 
errors are minimized in Sees. IV and V. The phase errors are minimized in Sees. VI and VII. 
The universal gates are tested numerically in Sec. VIII by the exact numerical modeling of 
quantum dynamics of the system. In Sec. IX we summarize our results. 



4 



II. FORMULATION OF THE PROBLEM 

The most important quantum algorithms exploit the property of quantum interference. 
In the result of implementation of these algorithms a small amount of "useful" states is 
amplified while all other states are suppressed. Hence, in order to implement a quantum 
logic a strict control over both moduli and phases of the amplitudes in superposition of states 
in the register of a quantum computer is necessary. For example, the quantum Control-Not 
(CN) gate CNj t fc transforms a superpositional wave function as 

2 L -1 2 L -1 

CN S al B Ml-i ■■■n i ...n k ...n Q ) = ^ B j\ n L-i ■ ■ ■ . . . n t © n k . . . n ), (1) 

3=0 j=0 

where © means sum modulo 2, n m = 0, 1, and L is the number of qubits. The CN gate for 
remote qubits can be implemented using a series of CN gates for the neighboring qubits, 
CN i+1 j or CNj (See Sec. |V111] of this paper or, for example, Ref. P7| .) Moreover, any 



quantum logic operation between remote qubits can be divided into elementary logic gates 
between neighboring qubits. Hence, in order to implement the universal quantum logic, it 
is sufficient to construct only the universal quantum gates between adjacent qubits. 
Usually, in real quantum systems, instead of Eq. ([[]) one has 

2 L -1 2 L -1 

CN ? ^2 B j\ n L-i ■ ■ - n i ■ ■ - n k ■ ■ ■ n ) = ^2 B 'j\ n L-i ■ ■ • n f . . .Ui®n k . . .n ), (2) 

3=0 3=0 

where the coefficients B'j in the right-hand side of Eq. (0) are slightly different from the 
coefficients Bj in the right-hand side of Eq. ([I]). This introduces the error in the quantum 
algorithm. A source of the error can be external, like noise, or internal, like unwanted 
transitions driven by pulses of a protocol. 

In this paper, we will consider only the latter errors. We will present protocols for 
implementation of basic quantum logic operations between the neighboring qubits which 
work for an arbitrary superposition of states in the Ising spin quantum computer. We use 
these operations to implement the CN gate between the end qubits, CN ,l-i of the spin 
chain, and estimate the phase and probability errors. 
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III. ISING SPIN QUANTUM COMPUTER 

The simplest Hamiltonian for the Ising spin chain placed in an external magnetic field 
can be represented as 

L-l L-2 „ L-l 

H n = -Y, " 2 J Wk+i ~ -y E { 4~ exp [-i + ^n)] + fc.c } = ^0 + V n (t) • (3) 

fc=0 k=0 k=0 

Here ft = 1, I£ is the operator of the z component of kth spin 1/2, 1^ = If ± 2^, cjfc is 
the Larmor frequency of the fcth spin, J is the interaction constant between the neighboring 
spins, Q n is the Rabi frequency (frequency of precession around the resonant transversal 
field in the rotating frame), v n is the frequency of the pulse, and cp n is the phase constant of 
the rath pulse, below called "phase" . The Hamiltonian (|3[) is written for one nth rectangular 
rf pulse. Below we omit the index n which indicates the pulse number. 

In the interaction representation, the solution of the Schrodinger equation can be written 
in the form 

*{t) = J2Cp{t)\p)exp{-iE p t), (4) 

p 

where E p and \p) are, respectively, the eigenvalues and eigenfunctions of the Hamiltonian 
H . 

Under the condition < J < u the pulse effectively affects only one kth spin in the 
chain Jl7j whose frequency Uk is close (near-resonant transition) or equal (resonant tran- 
sition) to the frequency of the pulse, v. In this approximation, the system of coupled 
differential equations for the coefficients C p (t) splits into 2 L ~ 1 independent groups. Each 
group consists of two equations of the form 

iC p = — (r2/2) exp[i(A pm t - <p)]C m , (5) 

iC m = -(fi/2) exp[-i(A pm t - (p))C p , 

where the states |ra) and \p) are related by the flip of kth spin, A pm = E p — E m — u, and 
we suppose that E p > E m . The solution of Eq. (|J) is 

C m (to + r) = [cos(A pm r/2) + i(A pm /X) sin(A pm r/2)] x exp[-?rA pm /2], (6) 

C p (t + t) = i(Vt/\ pm ) sin(A pm r/2) x exp[it A pm + z(rA pm /2 - (p)). 
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Here t is the time of the beginning of the pulse, A pm = A 2 pm + Q' 2 , r is the duration of 
the pulse, and the initial conditions are 

C m {t ) = 1, C p (t ) = 0. 

The solution for the transition from the upper state to the lower state is 

C m (t + r) = i(tt/\ pm ) sin(A pm r/2) x exp[-it A pm - i(r A pm /2 - cp)), 



C p (t + r) = [cos(A pm r/2) - i(A pm /\) sm(A pm r/2)] x exp[irA pm /2], (7) 

C m {t ) = 0, C p {t ) = 1. 

In the table below we present all states and their energies for the spin chain of three qubits 
with the Larmor frequencies uj — w, U\ — w+5u), u 2 = w+25uj (w is the Larmor frequency of 
the zeroth qubit) and with the Ising interaction constant between the neighboring qubits J, 



|0> = |0 2 0i0 >, E Q = -\w - \8u - J, |1) = |001), Ex = -\w-\5u, 

|2) = |010), E 2 = —\w — \5u + J, |3) = |011), Ez = \w-\5u, 

|4) = |100), E 4 = -\w + \5uj |5) = |101), E$ = \w+ \8u + J, 

|6) = |110), E & = \w + \5uj |7) = |111), E 7 = \w + \5uo- J. 

We will use this table for illustrative examples in the following sections. 
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IV. THE 2vr£;-METHOD 



Suppose, for example, that before the action of the pulse we have the superposition of 
states |0) and |4) in Eq. (f|). In order to organize, for example, the transition |4) — > |6), 
associated with the flip of the first (middle) spin, we apply the 7r pulse [Qr = it in Eq. 
(P)] with the resonant frequency v — E e — E± — w + 8u. We also have the near-resonant 
transition (or transition with small detuning from the exact resonance condition) |0) — > |2) 
with the detuning A2,o = E 2 — Eq — v = 2J which creates the unwanted state with the 
probability [see Eq. @] 



2,0 



.^2,0 



sin 



2Q 



(9) 
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where A2,o = y Al + ^ 2 = V4J 2 + Q 2 . In order to suppress this transition we take the 
Rabi frequency in the form 

2 / 



(See the 2irk method in Ref. |12|].) Then the argument of sinus in Eq. @ turns to zero, 



and the unwanted near- resonant transition |0) — > |2) is completely suppressed. 

V. GENERALIZED 2nk METHOD 

Next, we will show that one cannot suppress all possible near- resonant transitions using 
only one pulse. Since detunings for different transitions are different, one 7r pulse in general 
case creates a relatively large error. In this section, we will demonstrate how to compen- 
sate this error using additional correcting pulse. This procedure can be characterized as a 
generalized 2nk method since this protocol suppresses all near-resonant transitions for all 
quantum states in a superposition. 

Suppose that we have three quantum states, |0), |4), and |5) in a register of a quantum 
computer. Using the resonant pulse we produce the resonant transition |0) — > |2) by flipping 
the first qubit. We also have two near-resonant transitions: |4) — > |6) with the detuning 
Ag,4 = —2J and |5) — > |7) with Aj^ = — 4J. In order to suppress both transitions the Rabi 
frequency Q should satisfy the system of two equations 



A 



6.1 



n 



VW 2 



A 



7,5 



n 



V4K 2 - 1, 



where k and K must be integer numbers. Dividing the second equation over the first one 
and taking into consideration that ^7^/^6,4 1 = 2 we obtain the relation K 2 = 4k 2 — 3/4. 
From this equation one can see that if k is integer then K cannot be an integer number. 
Hence, if we suppress, for example, the transition |4) — > |6), then we get an error for the 
transition |5) — > |7), and vice versa. 

To proceed, we introduce the following notations. The pulse P® indicates the n pulse 
with the frequency z/ 4 00 which is resonant for the transition associated with the flip of ith spin 
when both its neighbors are in the states |0). In a similar way, we define the pulses P/° and 
P} 1 . The pulse P°° will be resonant, for example, for the transitions |0) — > |2) and |2) — > |0) 
[see Eq. (§)]. We define T®° as the transition (not necessarily resonant) associated with the 
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flip of ith spin when both its neighbors are in the states |0). For example, the pulse P™ 
gives rise to the resonant transitions T®° and two types of the near-resonant transitions, 
and T/ 1 . In general case of an arbitrary number of qubits, one notation, for example, T®°, 
means the transitions for all quantum states with the identical groups | . . . 0j + i0j0j_i . . . ) 
and | . . . 0j + ilj0j_i . . . ). So, now we discuss the types of transitions, but not the transitions 
between definite states. We should note that u} = z/? 1 , so that we use the same notation, 
P. 10 , for the pulses P/° and PP 1 . 

Now we continue the discussion on the 2irk method. All near-resonant transitions Tf° 
and T/ 1 , caused by the pulse P/ , can be suppressed using only one pulse, since the moduli 
of detunings are the same and equal to A = 2J . In this case, we take 

A = vw^i, (ii) 

where k is the integer number. 

As was shown above, pulses Pf° and P/ 1 generate errors because the detunings for tran- 
sitions Tj 00 , T/ and T} x are different. To correct this error, we suppress the unwanted 
transition T/ 1 for the pulse P®° and the unwanted transition T®° for the pulse P/ 1 with 
the moduli of detunings A 2 = 2A = 4J by the 2irk method with the Rabi frequency f2 2 
satisfying the expression 

= sjAkl - 1. (12) 

Below we take fc 2 — k and f2 2 = 2fi, Qi = Q, where Q is a parameter. Thus, we choose the 
value of the Rabi frequency for the pulses P/ 1 and Pf twice larger than the value of the 
Rabi frequency for the pulse P/ . 

The pulses Pf° and P/ 1 with the parameters (|T2"D generate errors in the result of the 
transitions T/ . In Appendix we demonstrate that the additional pulse with the frequency 
v}° can be used to correct this error. The real and imaginary parts of the amplitude of the 
created unwanted state are corrected by the proper choice of the duration and the phase 
of the correcting pulse. The correcting pulse does not generate errors since the moduli of 
detunings for unwanted transitions, Tf° and T/ 1 , are the same, so that these transitions can 
be suppressed simultaneously. 
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VI. UNWANTED PHASES 

In spite of the fact that all pulses described above are probability-corrected, they generate 
unwanted phases which should also be compensated by the protocol. In this section, we 
calculate all unwanted phases which appear in the quantum computer in the result of action 
of the rf pulses. 

Even in the case when all unwanted near-resonant transitions are suppressed, the un- 
wanted phases appear because for definite transitions A pm ^ in Eqs. (^) and (|7|). The 
problem to compensate those phases becomes complicated, since each pulse generates dif- 
ferent phases for different states. In this section we will show that the unwanted phases 
can be compensated by choosing the proper phases of the pulses of the protocol. Each logic 
operation requires its own set of phases of the pulses. However, if one has a set of phase- 
compensated universal gates, introduced below, one can realize a quantum logic using these 
gates as building blocks. 

Further analysis is conducted in terms of the probability-corrected pulses, Q™" n (ip), where 
m, n = 0, 1, i = 0, 1, . . . L — 1, and (p is the phase of the pulse. The pulse Q}°(<p) indicates 
one 7r pulse. It coincides with the pulse P/°, introduced above. The combined probability- 
corrected pulse, Qi°((p) [or consists of one n pulse Ff°, (Pi 1 ), and one correcting 
pulse. The correcting pulse does not implement quantum logic, it only removes unwanted 
states from the register of quantum computer. 

In the Table I we present the phases (see Appendix ||) acquired by different states (which 
can be in a superposition in a register of quantum computer) generated in the result of the 
different kinds of the probability-corrected pulses, Q™ n (ip). 

VII. PHASE AND PROBABILITY-CORRECTED UNIVERSAL GATES 

In the Ising spin quantum computer the parameters which allow one to compensate the 
unwanted phases are the phases of the rf pulses. From Eq. (Q) one can see that the phase 
(p of the pulse changes the phase of the wave function. In this section we provide the phase 
and probability correct protocols for realization of universal quantum gates on arbitrary 
superposition of states in the Ising spin scalable quantum computer. 

Our task is to construct phase-corrected quantum gates from the elementary probability- 
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state 


acquired phase 








| • • • Oi+iOiOi- 


-1 • • • ) 


-9 


tt/2 - v? + 7* 


-0-7 


... Oj+iljOj- 


-1 ■ • • ) 


e 


tt/2 + Lp - 7* 


+ 7 


| . . . lj+iOjOj- 


-1 • • • / 


tt/2 - <p* 


7T + 0/2 + 9 


7T-0/2-9 


• ■ • Ij+lliOj- 


_!...) 


Vr/2 + if* 

I 1 T 


7r-6>/2-e 


7T + 0/2 + 9 


... Oj + iOjlj_ 


-1 • • • ) 


71-/2 - Lp* 


7T + 0/2 + 9 


7T-0/2-9 


... Oj + iljlj_ 


_!...) 


tt/2 + Lp* 


tt-0/2-9 


7T + 0/2 + 9 


... Ij + i0jlj_ 


-1 • • • ) 


e 


+ 7 


7r/2 — Lp — 7* 


• • • Ij+llilj- 


-1 ■ ■ ■ ) 


-9 


-0-7 


tt/2 + Lp + 7* 



TABLE I: Phases generated by the probability-corrected pulses, QJ 171 ^). The asterisk indicates 
that the resonant transition from the state shown in the first column of the table to the other state, 
associated with the flip of the ith qubit, takes place. The phases, 0, 9 and 7, are defined in the 
Appendices [A| and |B|. 



corrected three-qubit elementary gates, Q™ 71 ^), m, n = 0, 1, i = 0, 1, . . . , L — 1. These gates 
act on the target qubit g^, for given configurations of its nearest neighbors, q^i and ft+i. We 
thus define the operations Q™ 71 ^) in terms of their effects on the states, | . . . g^igi^-i . . . ) 
as 

. for Qf{ip)=Qf{ip), 

I ... i+ i 9< 0i-i . . . ) - exp (-0)]|O m gA-i • • • ), 

I ••• i+ ifcli_i ...)-► exp i [tt/2 + (-1) 9 * (-if)} I ... (Vnftli-! . . . ), 

I ••• h+iQiOi-i ...)-► expi [tt/2 + (-p)] | . . . I m g70^i . . . ), 

I ••• li+iftli-i ■••>-> expi[(-l)« (0)]| . . . 1,+ig^-i • • • ), (13) 



for Q°°(^) 



Oi+iftOi-i . . . ) -> expi[7r/2 + (-l) ft (-^ + 7 )] I ...0 <+ i50i_i...), 
Oi+ifili-x . . . ) -> expz [tt + (— (0/2 + 6)] I . . . O^li-i), 
U+iqiOi-i ■ ■ ■ > -> exp? [tt + (— (0/2 + 6)] | . . . li+ifcOi-i . . .), 
U+^ih-! . . . > -+ expi [(-1) 9 < (0 + 7)] I ••• O^i^O^i . . . ), (14) 
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for Q} 1 ^), 



I • • • Oj+i^Oj^i 
I • • • Oi + igili_i 

| • • • U+lQi^i-l 
I • • • U+lQih-l 

One can see that each Qf' n pu 



expi[(-l)*(-0- 7 )]|...O i+1 g i O i -i...)> 
exp i [tt + (-I) 9 ' (-0/2 -6)]|... O i+ iftli-i • 
expi [tt + (-l) q > (-9/2 - 0)] | . . . li+iftOi-i . 
exp i [it/2 + (—l) 9i {-<p - 7)] I • • • h+iWih-i ■ ■ 



(15) 



lse introduces three different kinds of phases into different 
states. These phases arise from different types of resonant and near-resonant transitions 
initiated by the pulses. 



A. Not gate 

Each pulse has one externally controllable phase, ip, which can influence different subsets 
of states. Using combinations of three pulses on qubit i, one can introduce at most three 
independent phases to correct different unwanted phases generated for different states in the 
register of a quantum computer. 

For the Not gate, the correct transformation for the amplitudes of the states can be 
implemented by the sequence Q\ 1 {y$)($ 1 {<P2)Q*i i {<pi)- (The order of implementation of the 
operators is from the right to the left.) The transformation which results from this sequence 
of operations is, 



• • • O m giO;-i • • • ) 

• • • 0j+i5jlj_i . . . ) 

• • • U+iQA-i ■ ■ ■) 
■ ■ ■ h+iQih-i • • • ) 



exp % [tt/2 + (-1) 91 + 2 7 + 26)} \... O i+1 m~i ■ ■ ■ 
exp % [n/2 + (— l) 9i (— y?2 + 6 + 20)] | ... O i+ iftli_i . . . ) ; 
exp* [tt/2 + (-lH-y?2 + 6 + 20)] | ... l t+1 m-i ■■■): 
exp i [it/2 + (-1)* (-y? 3 + 26)} | . . . l^U-i ...). 



(16) 



It is easy to see, that the implementation of the Not gate requires the following set of 
phases: <fi — 27 + 29, y?2 = 6 + 20 and ip 3 = 29. The overall phase factor is equal to tt/2. 



B. Control-Not gate 



For qubits with homogeneous coupling between them, the correct CN a b gate, where a is 
the number of the control qubit and b = a ± 1 is the number of the target qubit, can be 
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pulse 


produced state 


acquired phase 


phase 


1 


Q n ( _ 5 _ 2t) 


... 0i + 2li+i0j0j^i . . . ) 


tt- ±0-9 


7T- ±0-9 


2 


Qi° (§0-e+ 7 ) 

T i V Z ' / 


| . . . 0j_|_2li+llj0j_i . . . ) 


±7T- (|0-9 + 7 ) 

A \ A 1 J 


|tt - 30 - 7 

Z ' 


3 


Q 4 10 (0) 


... 0j+2lj+i0j0j_i . . . ) 


¥ 

A 


-30-7 


4 


Q°° (|vr + 20- 
46 + 2 7 ) 


... 0j_|_20j_|_i0j0j_i . . . ) 


ivr+lvr + 20- 

Z 4 

49 + 2 7 


|tt - - 49 


5 


Qi+i (!*•) 


| . . . 0i+20j-i_i0j0j_i • . . ) 


-0 


|tt - 20 - 49 


6 


" t~r -L V 4 / 


| . . . 0j_(_20j_|_i0j0j_i . . . ) 


-0- 7 


|tt - 30 - 49 - 7 

4 ' 


7 


Q°° (-29) 


| . . . 0j+20i+ili0i-i . . . ) 


i^ + 29 + 7 

Z ' 


|vr - 30 - 29 


8 


Qf (-§0 + 9-7) 


... 0j + 20j + ili0j_i . . . ) 





In -20- 29 

4 


9 


Q 4 10 (o) 


| . . . 0j-(_20j+i lj0j_i . . . ) 





|vr - - 29 


10 


Q™! (20-49 + 2 7 ) 


... 0j + 20j + ilj0j_i . . . ) 


7T+ ±0 + 9 


fvr-±0-9 


11 


Qllx (0) 


... 0j + 2li+llj0j_i . . . ) 




fvr-±0-9 


12 


(o) 


... 0j + 2lj+lli0i_i . . . ) 


7T+ ±0 + 9 


It 



TABLE II: Modification of the phase of the initial state | . . . 0i+2lj+i0j0j_i . . . ) in the result of the 
action of the CNj+i^ gate pulses. 

implemented by the following sequence, 

CN a , b = Ql\0)Ql\0)QT(V8)Q b \0)Q b \^)Q™(^) 

^ 1 (^)g^ (^4)g° (^)^ 1 (o)g^ 1 (^2)^ 1 (^), (no 

where the phases (p n , are computed in Appendix [C], 

(p x = -56 - 2 7 , <^ 2 = §0 - + 7, y> 3 = f vr + 20 - 46 + 2 7 , ip A = fvr, 

^5 = fvr, ^ 6 = -26, ^ = -§0 + 0-7, y? 8 = 2# - 46 + 2 7 . 

In order to illustrate the action of the CN gate, in the Table II we show how each pulse of 
the gate CNj +l i modifies the phase (initially equal to zero) of the state | . . . 0j +2 li+i0i0j_i . . . ) 
due to the Table I. In a similar way, one can show that other states acquire the same phase 
7r/4 after the action of the CN gate, while their phases in the middle of the CN gate protocol 
can be different. 
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C. Not and Control-Not gate on edge qubits 

In order to make our set of elementary quantum logic operations complete, we also present 
the protocols for the logic gates on the edge qubits. The Not gate is implemented by the 
sequence of two pulses: 

Not, = Q}(0)Q?(0), (18) 

where i — 0, L — 1. In Eq. (|T8l) we suppose that the notation Q™({p), m — 0, 1, with one 
upper index means one ir pulse with the corresponding resonant frequency and the phase 
<f. The overall phase for the gate is ir/2. The gate CN a b with the edge target qubit b 
(b = or L - 1) is 

cN a , fe = Qi'mi'mTmi 1 (\A qi° (l^ Q° a ° Qtmn-mK-m, m 

with the overall phase — ?. The gate CN ai b, with the edge control qubit a (a = or L — 1), 
is 

CN a , b = Ql(0)Q°(0)Q 6 10 (0)Q 6 10 (0)Q°° (-60 + 26 - 2 7 ) Q l a + ^ - + 7 ) x (20) 

Q° a (f* - ^ + - 7 ) Ql°(0)Ql°(56 - 20 + 2 7 )Q 6 11 (-20) , 
with the overall phase f . 



D. A single qubit rotations around x axis for an arbitrary angle 

We denote the pulses for rotation of ith spin around x axis for the angle pn/2 as Q^i^p), 
m, n — 0, 1, i = 1, . . . , L — 2 for intermediate qubits, and Q™((p), i = 0, L — 1 for the edge 
qubits. Note that the case p = 1 was analyzed before. From the 27rfc-condition for Q p one 
obtains 

p = iry/k*-p?/4, «p = fv^ 2 + 3p74, (21) 

6> 7T 

tan p = — 2— tana p , tan/3 p = tana p cos0 p (22) 

oi p 2 gl p 

and 

7p = v /( 7r A;)2-( 7r + / 5 p )2 ) (23) 

where we took the same integer value, k p = k, as before. The pulse Q} p (f) is the single 
pulse with the Rabi frequency f2i iP = Ap/ 4k 2 — p 2 = f2 p , duration r p = pir/{l p , frequency 
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v 



10 



and phase p>. Note, that the Rabi frequency, Q p , decreases with p decreasing. The 



pulse Q}p(<p) [or Q®p((p)] is composed of the pulse with the Rabi frequency fl 2 , P = 2f2 
duration r = p7r/f2 2 ,p, frequency v\ x (^ 00 ), and the phase p, and a correcting pulse with the 
Rabi frequency Q CjP = A( / 9 P + ir)/jp, duration r CjP = 27 P /A, frequency Z7 10 , and the phase 
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—Op + p — At — 6 P (Pe p = @p + V 9 + Ato + The phases generated by the pulses 
for Q^if) are the same as the ones indicated in the Table I, if we replace 9 with 6 p , with 
G p , and 7 with j p . The resonant transitions for an arbitrary angle create a superposition 



of states. In the Table [TTj we specify the phases of both, excited and initial states after the 
transition. 



Pulse 



Qlfa) 

Qlfa) 
Qlfa) 
Qlp(<p) 
Ql x M 



Initial state 



0j+i0j0j_i . . . 
Oi+iljOj-i . . . 
0j+i0jlj_i . . . 
0j+iljlj_i . . . 
l«+i0j0j_i . . . 
Ij+iOjlj-i . . . 
Ij+iOjlj-i . . . 

li+llilj-l • • • 



Phase of initial state 



n/2 + p- 7 P 


n/2 + ip 


n/2 + p 
7p 

7r/2 + ip + 7 p 



Phase of excited state 



7r/2-p + 7p 

Ip 

it/2 - p 


it/2 - p 


n/2-p- 7p 



TABLE III: Phases acquired after a Qf™ pulse for resonant transitions 



The probability and phase correct rotation of a single qubit can be implemented by the 
operators: 



(24) 



Ufa) I . . . 0, . . . ) = cos(p7r/2) |...0 i ...) + ie ilp sin(p7r/2) | . . . 1, . . . ), 
Ufa)\ . . . Ij . . . ) = cos(pvr/2)| + ie~ itp sin(p7r/2)| . . . 0j . . . ). 

The operators Uj(cp) can be build using the following sequence of pulses Q™ n (<p): 
for intermediate qubits: 

U !^P) = GJ»9S°(2(7 + 26 + p + 1 p + e p ))9S°(0)gg(-47 - 86 + p - 2 lp - 26 p ) 

x Qf (-2(7 + 29 + 7p + ^p))gf (0)Q J 1 p (4^ + ^)Qj (0)gf (0); (25) 
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for the edge qubits: 

u;{<p) = Q]»Q°> + 2e p )Q)(-e p )Q 3 (d p )Q)(o)Q° 3 (o), (26) 

with the overall phase tt. 

VIII. NUMERICAL RESULTS 

In the above sections we were concerned only with the errors generated in the result of 
the near-resonant transitions with frequencies close to the frequency of the external field, 
when E p — E m — v ~ J. These transitions are associated with the flip of the resonant kth 
spin whose the Larmor frequency, Uf., is close to the frequency of the external field, LOf. ~ v. 
In general case, the near-resonant transitions generate large errors. In the preceding sections 
the procedure is described how to suppress all these transitions. However, since the external 
rf field affects all spins in the system, there are non-resonant transitions associated with 
flips of other non- resonant k'th spins. Because the non-resonant transitions \p') — > \m') have 
large detuning, E p i — E m > — v ~ 5u\k — k'\, their probabilities are small and proportional to 



fj?/\h - k'\ 2 , where [0, [K 

fi = (n/26u). (27) 

For typical parameters fi ~ 10~ 4 . Since in the system there are only non-resonant transitions 
(and all near- resonant transitions are suppressed), for a definite protocol fi is the only 
parameter which defines the errors in our quantum computer. This is demonstrated in the 
numerical simulations below. 

In order to demonstrate the action of the logic gates on superpositional states, described 
in this paper, we simulated the exact quantum dynamics [l?]] of the system with the Hamil- 



tonian (|]) during implementation of the CN gate (0) between the edge qubits of the spin 
chain, CN ,l-i, where the Oth qubit is the control qubit and the (L — l)th qubit is the target 
qubit. (In a similar way one can implement any other CN gate CNjj, where i, j — 0, . . . , L— 1 
and i ^ j, for example CN L _ 1>0 .) Using the Swap operations, 

Si,i+i = Sj+i^ = CNj j j +1 CNj +lj jCNj i j +1 = CNj +lj jCNj ii+1 CNj +li j, (28) 

we move the control Oth qubit to the L — 2th position, implement the CN gate CN l _ 2> l-1) 
and using the same Swap gates we return the control qubit to its initial Oth position. The 
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FIG. 1: Deviation of phases from the common phase for different states of the superposition after 
implementation of the CNo,z,-i gate, L = 7, f2 = flk [see Eq. (|10|)], k = 2, J = 1; (a) 5u> = 10 4 , 
(b) 5uj = 2 x 10 4 . 



whole procedure can be written in the form (read from the right to the left) 

CNq,£_1 = OOjlSi^ • • • Sl_4 i l_3Sl_3 j l_2CNz / _2 i L_iSl_2,L-3Sl-3,L-4 • • • 82,181,0- 



(29) 



For simulations we initialized the system in the superpositional state (P with randomly 
chosen real positive values of the coefficients Bj, subject only to the normalization condition 
Yl"j=o 1 Bj ~ 1- Since all Bj are real and positive the phases of all states of the superposition 
are the same and equal to zero. (We should note that our protocols work for any set of 
initial phases.) 

In Figs. l(a,b) we show the deviation of phases for different states, (pj— <&, j — 0, . . . , 2 L — 1 
[here $ is the total phase and tpj = arctan[Im(i? / :; )/Re(i? / J '))], where Im and Re stand for 
the imaginary and real parts of the amplitude B'j in Eq. (§)], from the total phase $ after 
implementation of the CNo,l-i gate for two values of Sou [and correspondent values of \i in 
Eq. (|27|) 1. From comparison of Fig. 1(a) with Fig. 1(b) one can see that decreasing twice 
the value of /i leads to decrease of the error in the phases of all states by the factor two. 
This supports the fact, mentioned in the beginning of this section, that for our protocols /x 
is the only parameter which defines the errors in the quantum computer. Physically, this 
means that only the non-resonant transitions contribute to the errors and all near-resonant 
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O 100 200 300 400 500 

Number of pulses 

FIG. 2: The phase error as a function of the number of pulses for four values of L. 5ui = 10 4 , other 
parameters are the same as in Figs. l(a,b). 

transitions are completely suppressed. 

In Fig. 2 we plot the phase error as a function of time. The phase error is defined as a 
the maximum max | — <&|, j — 0, . . . , 2 — 1. In Figs. l(a,b) the phase error is equal to 
the distance between the two dashed lines. During implementation of each elementary CN 
gate between the neighboring qubits (using 12 pulses for intermediate qubits) the phases 
of different states of the superposition are different, but after implementation of each gate 
CNj ji+ i or CNj+i^ the phases are equalized, since our CN gates between the neighboring 
qubits are probability and phase-corrected. That is why in Fig. 2 the phase error is plotted 
with the interval of 12 pulses for intermediate qubits and 9 or 10 pulses for the edge qubits. 
As follows from the results of the calculations, the phase error is mostly generated by the 
pulses Q®° and Qj 1 . From Fig. 2 one can see that the phase error grows linearly with the 
number of these pulses. For larger number of qubits in the chain the error is larger, because 
the implementation of the protocol requires larger number of pulses. [The total number of 
pulses Q™ n is 2 x 36 x (L - 2) - 5.] 

The probability errors for the gate CN 0j l_i are defined as the moduli of the difference 
Pj{L) = \ \Bj\ 2 — \B'j\ 2 \, where the coefficients Bj and B'j in Eq. (0) are computed for the 
systems of L qubits. The values of Pj(L) are shown in in Fig. 3, and the relative errors 
Pj(L)/\Bj\ 2 , are plotted in Fig. 4 for different states of the superposition. From Fig. 3 
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FIG. 3: The probability of errors Pj(L) for different states of the superposition and for different 
values of L. 5lo = 10 4 , Vt = Vt k [see Eq. (0)], k = 2, J = 1. 

one can see that the absolute values of probability errors Pj{L) decrease with L increasing, 
since the initial occupation probabilities are smaller for larger values of L (because of the 
normalization condition Ylj^ 1 \Bj\ 2 = 1)- On the other hand, as Fig. 4 demonstrates, the 
relative probabilities of error increase with L increasing because of increasing the number of 
pulses in the protocol. 



IX. CONCLUSION 

In this paper the first attempt is made to develop the protocols for implementation of 
the universal quantum gates on an arbitrary superposition of quantum states in a multi-spin 
quantum computer. The quantum gates are constructed for a homogeneous spin chain (the 
spins and interactions between them are identical). In spite of the fact that such a computer 
is not realized experimentally (because of severe technological challenges), it provides a 
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FIG. 4: The relative probability of errors Pj(L)/\Bj\ 2 for different states of the superposition and 
for different values of L. Other parameters are the same as in Fig. 3. 



relatively simple model with basic features common for many solid-state scalable quantum 
computer proposals. 

We solved the general problems which one should face working with superpositional states 
in different types of quantum computers characterized by the constant (not switchable) inter- 
action between qubits (for example in Kane's quantum computer fl3fl). Since the transition 
frequencies required to flip the jth qubit in the states | . . . Oj + ig,,0.,_i . . . ), | . . . Ij + igj0j_i . . . ) 
(or | . . . Oj +1 qjlj_i ...)), and | . . . lj + iqjlj^i . . . ) (qj = 0, 1) are different, one inevitably will 
create an error by nipping the jth qubit. For an Ising spin quantum computer we corrected 
this error using additional correcting pulse which removes unwanted states from the register 
of quantum computer, created in the result of the unwanted transitions. 

Another important problem solved in this paper is a minimization of the phase errors. 
Even in the situation, when the probability errors are minimized, each pulse of the protocol 
produces different phases for different quantum states of the superposition. In the Ising spin 
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quantum computer the phases can be equalized by proper choice of phases of the rf pulses. 
Each protocol (for CN gate, one-qubit rotation, and others) has its own set of phases. 

The phase and probability errors studied in this paper, decrease when the Rabi frequency, 
Q, decreases or the frequency difference between the neighboring qubits, Sou, increases. In 
particular, all errors tend to zero when Q — > or 5uj — > oo. 

Our simulations show that the phase error less than 0.1 radians for the studied protocols 
can be achieved for the ratio 8u/J> 10 4 . In our example with phosphorus impurity donors 
in silicon (see Introduction) the ratio 5u/J is about 170. We will mention three ways to 
increase this ratio. The first way is to increase 5u using even greater gradients of the 
magnetic field suggested in |/J. The second way is to decrease the interaction constant J. 
This can be easily achieved by either increase of the distance between the impurity atoms 
or using the atoms with smaller effective radius. However decrease of J will reduce the 
clock speed of a quantum computer below MHz region. The third way is to develop more 
sophisticated sequences of electromagnetic pulses, which further suppress the phase error 
caused by non-resonant transitions. 
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APPENDIX A: CORRECTING PULSE 

Suppose that only the state \m) — \ . . . Ij + i0j0j_i) is initially populated, C m (to) = 1. We 
apply the pulse P} 1 with the frequency v\ x = uoi — 2J, where u>i is the Larmor frequency 
of the ith qubit. The detuning for the transition |m) — ► \p) = | . . . lj + il,;0j_i) is A pm = 
E p — E m — v\ x = A, where A = 2.1 . From Eq. @ after the pulse with the frequency v\ x 
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and phase $ we have 

C m = [cos a(k) + if(k) sin a{k)] e - im/2 (Al) 

C p = ig(k) sina(A;)e i[e(fc)/2 - ,I>+Ato1 . 

The Rabi frequency Q 2 of this pulse is chosen to suppress the transitions T®°. In this case 
the quantities a, 9, f, and g, defined below, are the functions of only one parameter k. From 
the 27r/c-condition one has 



^ = sjAkl - 1, A 2 = 2A. (A2) 

We choose 

Q 2 = 2Q, k 2 = k. (A3) 
For the transition T 10 one has A pm = A in Eq. ( |A1|) , so that 



a «=?V^=in/(£) H-l = iVFT374, (A4) 
where r 2 = 7r/f2 2 is the duration of the first pulse. The other quantities are 



A 



9{k) = Ar 2 = ^ = JV^l, (A5) 
ii 2 ^ 



/w = vWW = V^' (A6) 

= ; ^ = ; 1 (A7) 

Va 2 + v^ 2 + 3/4 

After the correcting pulse with the frequency v 10 , phase tp 11 (here the superscript indicates 
that the first pulse has the form Pi 1 ), and the duration r c , one has 

C m (t) = C m cos + iC p sin pe iip '\ (A8) 

C p (t) = C p cos p + iC m sin pe~ iip1 ^ , 

where 

t = t + t c , p= (A9) 
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From the conditions Re{C p (t)} = and lm{C p (t)} = 0, where Re and Im stand for, respec- 
tively, the real and imaginary parts of C p (t), one obtains the following equations: 

tan0 = — /tana, y^ 1 = — 9 + $ — Ato — ©, tan/? = — g tan a cos 0, (A10) 

where we do not indicate dependence of parameters on k. The second and third equations 
in ( |A10| ) define, respectively, the phase ip]} and the duration of the correcting pulse required 
to correct the error. 

If one considers the transitions generated by the pulse P} 1 in the opposite direction 
| . . . Ij + ilj0j_i . . . ) — > | . . . lj+iOjOj-i . . . ), one should just change the signs as [see Eqs. (||) 
and (0)] 

/->-/ 0->-0, Ato^-Ato, <f c ^-<Pc (All) 

in Eq. QA10 ). After this transformation both sides of the second Eq. ( |A10|) change their sign, 
so that the expressions for the parameters of the correcting pulse, tpl 1 and (3, do not change. 
Hence, one can suppress both unwanted transitions, | . . . I i+ i0j0j_i ...)—> | . . . I i+ ilj0j_i . . . ) 
and | . . . Ij + ilj0j_i ...)—> | . . . Ij + i0j0j„i . . . ), simultaneously. 

For the transitions generated by the pulse Pf° one should change the sign of A, so that 
the expression for the phase ip® of the correcting pulse takes the form 

<p™ = e + $ + At + e, (Ai2) 

where is defined in the first equation ( |A10| ). 

For k > 1 the value of (3 is small and negative. We take the duration of the pulse 
P* = (3 + 7T, so that C m {t) in Eq. (|A8[ ) changes its sign, and C p (t) is again equal to zero. 
After the correcting pulse one has 

C m (t) = - exp(-i0/2 - i6), C p (t) = 0. (A13) 

One can suppress the non-resonant transitions, T®° and T} 1 generated by the correcting 
pulse with moduli of detunings A, choosing the value of Q c satisfying 27r£;-condition, 

Xr ' nk c , A c = VA 2 + ^, r c = ^, (A14) 



2 ' v ' O 



where fc c is the integer number. From Eq. ( |A14j) one has 



A Ifirk,^ 1 



a c nrJ~ L (A15) 

Below we take k c = k. Since (3* ~ n, one has Q c ~ 2fii « f^2 where Qi and are defined 
in Eqs. flll]) and (0). 
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APPENDIX B: UNWANTED PHASES 



We calculate here unwanted phases generated by the probability-corrected pulses Q^itp), 
Q®°(<p), and Qj 1 ^). From the first equation in (|A13|) the pulse Q} 1 generates the phase 
n - 9/2 - Q for the states | . . . Ii + i0j0;-i . . . ) and | . . . O^iO^-i . . . ). From Eq. QATTD 
the same pulse generates the phases n + 9/2 + 6 for the states | . . . Ij + ilj0j_i . . . ) and 
| . . .Oj + iljlj_i . . . ). In a similar way one can deduce the phases generated by the pulse 
Q®°(ip) for the states | . . . l i+ i^0^_i . . . ) and | . . . Oj + igjlj_x . . . ), where = 0, 1. 

Now, based on the results of Appendix we calculate the phase generated by 

the pulse Q} 1 ^) for the state | . . . Oj+iOjOj-i) Since detuning for the transition 

| . . . 0j + i0j0j_i ...)—> | .. . 0i + ilj0i_i . . . ) is positive and equal to A 2 = 2A, from Eq. (H) the 
first 7i pulse [of the combined pulse Q} 1 ^)] generates the phase — A 2 /^2 = — A/fi = —9. 
The duration of the correcting pulse with the frequency is defined by the third equa- 
tion in Eqs. ( |A14j) , and the detuning is A. From Eq. (|6|) the phase acquired by the state 
| . . . 0j + i0j0j_i) ... in the result of action of the correcting pulse is 

A ' n -/5*^ = -v / (^ c ) 2 -/3* 2 = -7, (Bl) 



2 n c 

where we used Eq. (|A15|) . After the action of two pulses, which constitute the probability 
corrected pulse Q} 1 ^), the phase acquired by the state | . . . 0j + i0j0i_i . . . ) is — 9 — 7. We 
should note that for both pulses (main and correcting) the transition | . . . 0j + i0j0j__i ...}—> 
I . . .0 i+ ilj0j_i . . . ) is suppressed by the 2nk condition. In a similar way one can deduce 
other phases for other pulses in the Table I. 



APPENDIX C: PROTOCOL FOR THE CN GATE 



Our aim is to implement the CN gate between neighboring qubits, CN a ^, using only 
operations which affect the control and the target qubits, q a and g^, respectively. Thus, for 
complete generality we need to consider the states of the control and target qubits, as well 
as their neighbors. The CN gate in a homogeneous spin chain without phase correction can 
be implemented by the following sequence 

Qa 1 (^)Qi Mgr(^)Q^ (^)Qi 1 (^)Qi (^3)g° (^)^ 1 (^i). (ci) 
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This sequence can be described as follows: (in terms of the transformations for the ampli- 
tudes, we note that the phases of the various states will not be equal) (i) a controlled flip of 
qb if both of its neighbors are in the "1" state, (ii) a Not operation on q a , (iii) a controlled 
flip of qb if both of it's neighbors are in the "0" state, and finally (iv) another Not operation 
on q a which returns it to its initial state. 

The transformation resulting from flCT| ) is (we indicate in this Appendix only four relevant 
qubits of the spin chain) , 



|0000) 
|0001) 
|0010) 
|0011) 
|0100) 
|0101) 
|0110) 

loin) 

1 1000) 
|1001) 
|1010) 
1 1011) 
|1100) 
1 1101) 
|1110) 

1 1111) 



:[-7-0/2+e-¥>2+^8]|QQQQ\ 
.[+7+6>/2-e- ( p 2 +<Ai] |0QQ]^ 

.[+j+e/2-e-ip3+<f 7 ] 1 00 10) 
\-i-d/2+e-ip3+v7] | ooii) 

;[7T/2+7+30/2-e+ V 3+^5-V6]|QllQ\ 

[7r/2+7+0/2+8+pi+p2-P6]|Q]^]\ 
:[7r/2-7-36»/2+e+v? 2 -<P5-V7] |0100) 
:[7r/2-7-e/2-e-^i+^3-^ 7 ] |Q ^Q]^ 

■.[—r-0/2+e-<p 3 +<p T ] 1 1000) , 

:[+7+0/2-0-v?3+(,£>7] I ]_QQ]^ 

:[+7+6»/2-e-^ 4 +¥'8] 1 1010) 

'■[-7-o/2+e-<p4+<p B ] I ion) 

:[7r/2-7-e/2+3e+ip4+V5-^7] |0H0) 
:[7r/2+7+e/2+e+ V i+^3-¥'7] 1 1 1 1 1) 

:[7r/2+7+e/2-3e+^3-V5-^s] I iioo) 

,[n/2-'y-0/2-e-tpi+tp4,-tp8] | HQ1) 



(C2) 



However, when this sequence of gates is used, it is impossible to make the phases on 
all of the states equal. This can easily be seen by examining, for example, the first two 
transformations in Eq. (|U2|) . Here, it's impossible to equalize the phases on the two final 
states, for arbitrary values of 9, 0, and 7. Thus, in order to produce a phase-corrected CN 
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gate, it's necessary to add extra pulses which can introduce more controllable phases ipj. 
A sequence that we have found to work is, 



Site), 



Ql l MQl°MQTM d(o)Q°LM QlLM, (C3) 



where the braces indicate additional operations not found in Eq. ( |C1|) . The transformation 
which results from this sequence is, 

|0000) -> e ^[T-7-50/2+e- V 3- w +^] |0000), 
|0001) -> e i[7r+7+50/2-8 ^ 2 ^ 3+V8] |OOOl), 
|0010) -> e ^[T+7+5e/2-8- V4 + w +H |0010), 

1 001 1) — > e *[ 7r -7-5e/2+e+v5 2 -¥'4+¥'9] |0011), 

1 100) — > e i l 37T / 2 + , y-o/' 2 + i P2+<Pi+v6-v>a] |ono) 

|0101) e i[3n/2+ r y+50/2+Q+<pi+'P3-<pr-<P8]\Qm\ 

|0110) — > e *[ 37r / 2 -')'+^/2+e-ip2+^3-^6-¥'9] i oxoo), 

|0111) — > e *[ 37r / 2 -7-56'/2-e-^i+ip4+V7-¥ 3 g] |0101) 

1 1000) -> e ik-7-e«/2+e-V4-9T+<P9] | iooo) , 

1 1001) — > e d 7r +f+56>/2— e-v?2— 1 1001), 

1 1010) — > e i [ 7r ^7-56 , /2+6-V5+V7+¥'io]|1010^ 

1 1011) — > e i[T+7+50/2-8 +V2 - V6 +^io] | lOH) 

1 1100) — > e *[ 37r /2-7-50/2+30+^ 2 +<^5+^6-¥'9] |0H0) 

| HOI) — > e *[ 37r /2+7+50/2+e+^i+ip4-¥'7-¥'g]|iiii) 

1 1110) — > e *[ 37i "/2+7+56'/2-30-v32+¥'4-V6-Vio] |H00) 

|1111) — > e ?; [ 37r /2-7-56'/2-e-(pi+ip5+¥'7-¥'io]|ll01) (C4) 



This transformation introduces more than enough independent variables, (pj, to equalize 
all phases. To find the correct values for the tpj's we set all the phases above equal to the 
same number, and solve the resulting system of equations. The solution in the form 
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ip x = -2 7 - 50, <p 2 = 7 + f - 6, <p 3 = f + 2 7 + 20 - 46 + 2y? 9 - Pio, 
V?4 = x + ¥5 = x + <pio, y?e = -26 + <^ 9 - Vio, ^7 = -7 - f + ©, 

V? 8 = 2 7 + 20 - 46 + 2^ 9 - <^ 10 , (C5) 

equalizes the phases of all states to 7r/4. There are only eight independent tp/s in these 
equations, so that (p 9 and ip w can be set to zero. Thus, the pulse sequence (|C3|) with the 
phases flCS]) implements a phase corrected CN gate. 
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